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1 Introduction 

The study of step 2 subRiemannian manifolds has the Heisenberg group as a prototype. This 
is a noncommutative Lie group with the base manifold R 3 and endowed with a nonintegrable 
distribution spanned by two of the noncommutative left invariant vector fields. This structure 
enjoys also the property of being a contact structure or a CR-manifold. The study of the 
subRiemannian geodesies on the Heisenberg group started with the work of Gaveau [9]. One 
trend in the literature is to use the geometry of the Heisenberg group to describe the Heisenberg 
Laplacian and its heat kernel, see Beals, Gaveau, Greiner [1,2,3,4]. Later, this structure led 
to generalizations of the Heisenberg group as can be seen in Calin, Chang, Greiner [5,6] and 
Chang, Markina [7]. For more fundamental issues on subRiemannian geometry, see Strichartz 
[10]. 

One question is, if there is another 3-dimensional noncommutative Lie group besides the 
Heisenberg group. The answer is positive and this paper deals with it. The sphere S 3 and the 
group SU(2) have these properties. More important, this is the first compact subRiemannian 
manifold has been under consideration. 

The present paper starts with an introduction to quaternions and then defines the 3- 
dimmensional sphere as the set of quaternions of length one. The quaternion group induces on 
S 3 a structure of noncommutative Lie group. This group is compact and the results obtained 
in this case are very different than those obtained in the case of the Heisenberg group, which 
is a noncompact Lie group. Like in the Heisenberg group case, we introduce a nonintegrable 
distribution on the sphere and a metric on it using two of the noncommutative left invariant 
vector fields. This way S 3 becomes a subRiemannian manifold. It is known that the group 
SU (2) is isomorphic with the sphere § 3 and represents an example of subRiemannian manifold 
where the elements are matrices. The main issue here is to study the connectivity by horizontal 
curves and its geodesies on this manifold. In this paper, we are using Lagrangian method to 
study the connectivity theorem on S 3 by horizontal curves with minimal arc-length. We show 
that for any two points in S 3 , there exists such a geodesic joining these two points. 
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2 The Quaternions Group Q 

William R. Hamilton introduced quaternions in 1873. They are obtained as a non-commutative 
extension of the complex numbers. Let Q be the set of matrices ( %; _ ] with complex number 
entries a, (3 € C. The set Q becomes a group with the usual matrix multiplication law, called 
the 9 r mp . The identity element is I . (J J) and the inverse of , . («, <L 

is q- 1 = \jj ~P) , where A = \a\ 2 + \f3\ 2 . Let 



Then 



i \ ._/0 l\ 

-J ' J ~l-1 ' [i o 



i 2 = j 2 = k 2 = —I, ij = -ji = k, 
jk=-kj = i, ki= ik = j. 



(2.1) 



Any real number aeR can be identified with ai by the one to one homomorphism ip : M — > Q, 
a > 
-a, 



a 

< < c ( a ) = ( n I . If a = ao + aii, /3 = &o + 6ii € C, then an elementary computation yields 



-o _ = a I + ail + 6 J + Oik. 

a) 

The set {I, i, j, k} is a basis of Q. Let q = ai + bi + cj + dk and q 1 = a'l + b'i + c'j + d'k be two 
quaternions. 

A straightforward computation shows 

qq' = (aa' - bb' - cc' -dd')l+ {ab' + ba' + cd! - dc')i 

+ (ac + ca + db' - bd')j + (ad' + da' + be - cb% 

The conjugate of q = ai + bi + cj + dk is q = al — bi — cj — dk. A simple computation shows 
that qq = a 2 + b 2 + c 2 + d 2 = A, q -1 = —q. The modulus of q is \q\ = \/A- Since i, j, k do not 
commute, the multiplication of quaternions is noncommutative. 
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3 § 3 as a noncommutative Lie group 



Consider § 3 = {q E Q; \q\ = 1}. For any q x ,q 2 G § 3 , \q x 1 q 2 \ = 
of Q. The law group on § 3 is 



= 1, i.e. S 3 is a subgroup 



(a, b, c, d) * (a', b' , c', rf') 
= (aa' - bb' - cc' - dd' , ab' + ba' + cd' - dc' , ac' + ca' + db' - bd', ad' + da' + be' - cb'), 

with the identity element (1, 0, 0, 0) and the inverse (a, b, c, d)^ 1 = (a, —b, —c, —d). In this way 
(§ 3 , *) becomes a compact, noncommutative Lie group. 

We shall work out the left invariant vector fields with respect to the left translation L q : 
S 3 — > § 3 , L q h — q* h. These are vector fields on S> 3 such that 

dL q (X h ) = X q « h , V<z,/ieS 3 . 

The set of all left invariant vector fields on S 3 

£(§ 3 ) = {X:dL q X = X,Vqe § 3 } 

is called the Lie algebra of § 3 . £(§ 3 ) is a subalgebra of C{Q), i.e. it is a subset closed under 
the Lie bracket [ , ]. 

In order to find a basis for £(§ 3 ), we shall find first a basis for £(Q). Let X be a left 
invariant vector field on Q. Then X q = (dL q )X e , where e = (1,0,0,0) is the identity of Q. If 
write 

4 q 
i=l 

the components are 

X\ = {dL q )X e { Xl ) = X e { Xl o L q ), i = 1, . . . ,4, 
where Xi is the i-th coordinate. 

Next we shall compute Xi o L q , i = 1,2, 3, 4. Let q — (a, b, c, d) and q' — (a', b', c', d'). Then 
{xi o L q ){q') = x^Lqq') = aa' - bb' - cc' - dd' . (3.2) 

Using 

a = x\(q), a' = x 1 (q'), b = x 2 (q), b' = x 2 (q'), 
c = x 3 (q), c' = x s (q'), d = x 4 (q), d' = Xi(q'), 
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equation (|3.2p becomes 

(x 1 o L q )(q') = x 1 (q)x 1 (q') - x 2 (q)x 2 (q') - x 3 {q)x 3 {q') - x 4 (q)x 4 (q'). 

Dropping the argument q' , yields 

xi o L q = Xi(q)xi - x 2 (q)x 2 - x 3 (q)x 3 - x 4 (q)x 4 . 

By a similar computation we find 

x 2 oL q = xi(q)x 2 + x 2 (q)xi + x 3 {q)x 4 - x 4 (q)x 3 , 
x 3 oL q = xi (q)x 3 + x 3 (q)xi + x 4 (q)x 2 - x 2 (q)x 4 , 
x 4 oL q = x\ (q)x 4 + x 4 (q)xi + x 2 {q)x 3 - x 3 (q)x 2 . 

The hrst component becomes 

X\ = X e {x x oL q ) 

= xi(q)X e (xi) - x 2 (q)X e (x 2 ) - x 3 (q)X e (x 3 ) - X 4 (q)X e (x 4 ) 
= Xl (q)Xl - x 2 (q)X 2 e - x 3 (q)X 3 e - x 4 {q)Xi. 

Dropping q and denoting = X l el i = 1,2, 3, 4, the previous relation yields 

x 1 = X1 e -x 2 e -x 3 f - X4 ti 4 . 

In a similar way we obtain 

X 2 = x^ + x^+x^-x^ 3 , 

x 3 = Xl e + X3 e + x£ - x 2 z\ 

X 4 = x^ + x 4 S} + x 2 £ 3 - x 3 £ 2 . 
Left invariant vector fields Z can be written as 

z =Y, zl ^- = {^e-x 2 e-x 3 e-x 4 e)-^ 

i—l 

+{x l e + x 2 e + x 3 e ~ x 4 e)^ 
+{x l e + x 3 e + X4 e ~ X2 e) ^ 
+( Xl e + X4 e + X2 e - X3 e) 
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Reordering, yields 



Z = 



x 2 



dx 2 



X 3 



e 
+e 
+e 



dx 3 



X4 



8xa. 



d 


d 


d 


d 1 


dx\ 


ox 2 


ox 3 


0x4 J 


d 


d 


d 


d 1 


X 3~X 

OX i 


~ x ^ — 

ox 2 


dx 3 


+ X 2 - 

OX4 J 


d 


d 


d 


d 1 


X4 "5 — 

OX\ 


+ ^3^ 

OX2 


dx 3 


0x4 i 



where 



N 



4 a 



is the normal vector field to § 3 , and 

X = x 2 

Y = x 4 



T = 



x 3 



d 


d 


d 


d 


dx\ 


~ Xl 7T~ 

0x2 


~ X4 a — 

ox 3 


+ x 3 - — 

0x4 


d 


d 


d 


d 


dxi 


~ X 3JT~ 

0x2 


ox 3 


~ Xl JT~ 
0x4 


d 


d 


d 


d 


dxi 


+ Xi ^~ 
0x2 


dx 3 


~ X2 jr~ 
0x4 



Using the 4 x 4-matrix representations: 



h 



1 

-10 



1 









1 














1 


-1 














-1 
















-1 



1 

-1 







(3.3) 



and the standard inner product ( , ) on K , we also can write 

N = (xU,V x ), X = -(xI u V x ), Y = -(xI 3 ,V x ), T=-(xI 2 ,V x ), 

where U is the 4x4 identity matrix. 

It is easy to see that {N, X, Y, T} is a basis for £(Q). In order to find a basis for £(§ 3 ) we 
shall consider the vector fields which are tangent to S 3 . Since a computation provides 

(N, X) = (N, Y) = (N, T) = 0, 

it follows that X,Y,T are tangent to S 3 . Hence any left invariant vector field on (§ 3 ,*) is a 
linear combination of X, Y and T. Since the matrix of coefficients 



x 2 


-Xi 


— X4 


x 3 


X4 


-x 3 


x 2 


-Xi 


x 3 


X4 


-X\ 


-x 2 
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has rank 3 at every point, it follows that X, Y and T are linear independent and hence form a 
basis of L(S 3 ). It worth nothing that the above rank property is not true without the constraint 

We also note that {X, Y, T, N} form an orthonormal system with respect to the usual inner 
product of K 4 . Then £(Q) = £(§ 3 ) © RN. 



3.1 The horizontal distribution 

Let H — span{X,Y} be the distribution generated by the vector fields X and Y. Since 
[Y,X] = 2T £ 7i it follows that H is not involutive. We can write £(§ 3 ) = H@WT. The 
distribution Ti. will be called the horizontal distribution. Any curve on the sphere which has 
the velocity vector contained in the distribution H will be called a horizontal curve. 

Remark 3.1 Notice that this situation differs from the Heisenberg group since the choice of 
the horizontal distribution is not unique. Because of [X,T] = —2Y and [Y,T] = 2X we could 
chose H = span{X,T} or H = span{Y,T}. The geometries defining by different horizontal 
distributions are cyclically symmetric, so we restrict out attention to the H — span{X, Y}. 

The following result deals with a characterization of horizontal curves. It will make sense 
for later reasons to rename the variables j/i = £3, j/2 = 2^4 • 

Proposition 3.2 Let 'y(s) = (xi(s),X2(s),yi(s),y2(s)) be a curve on § 3 . The curve 7 is 
horizontal if and only if 

(x,y) = (x,y). 

Proof: Since {X, Y, T} span the tangent space of the sphere S 3 we have 

7 = aX + bY + cT. 
Using that {X, Y, T} is an orthonormal system on § 3 , we have 

c = (7,r> 

= ixyi + x 2 y 2 - yixi - 1)2X2 
= (x,y)-(x,y). 

The curve 7 is horizontal if and only if c = i.e. (x, y) = (x, y). 

■ 

Hence a horizontal curve 7 = (x\, X2, yi, 1/2) : (0, r) — > R 3 has the velocity 7 = aX + bY 
where 



a 



b = 



(7, X) = ±ix 2 - x 2 xi - yiy 2 + myi 
(7, Y) = iiy 2 - x 2 yi + yix 2 - ij2X\. 



(3.4) 
(3.5) 



G 



If we consider a metric on H such that {Y, X} are orthonormal vector fields, then the length 
of the curve 7 is 



£(7)= f ^a(s)2 + b(s)2ds. 
Jo 



The horizontal curves minimizing length will be treated in a next section. 

In the rest of this section we shall make a few considerations regarding the horizontal 
distribution given as the kernel of a one-form. 

Consider the following one-form on R 4 

u> = Xidyi — y\dxi + X2dy 2 — yidx2 = xdy — ydx. (3-6) 

One can easily check that 

u(X) = -x 1 y 2 -yiX2+x 2 yi+x 1 y2 = 
u(Y) = xix 2 - y\yi - xix 2 + y\yi = 
w(T) - -x\- y\-x\-yl = -1^0 
lo(N) = xiyx - yixi + x 2 y 2 - V2X 2 = 0, 

where 

N = x 1 d Xl + x 2 d X2 + yid Vl + y 2 d V2 

is the unit normal to the sphere § 3 . Hence kerw = span{X, Y, N} and the horizontal distribu- 
tion can be written as 

S 3 3 p -> 7i p = kevuj n T p § 3 . 
Hence a vector v = (v, w) — (vi,V2,wi,w 2 ) € H if and only if 

(x,v) + (y,w) = 0(-<=>v-LN) 
(v,y)-(x,w) = = 0) 

\x\ 2 + \y\ 2 = lM(x,!/)e§ 3 ). 



4 Connectivity theorem on S 3 

The present section deals with the global connectivity property by horizontal curves. Even if a 
result of Chow [8] states the connectivity by piece- wise horizontal curves, here we shall prove a 
globally smooth version of it. More precisely, the main result of this section states that given 
two points on § 3 , there is a horizontal smooth curve joining them. 

Before proceeding to the proof, we need to introduce some terminology and notations. This 
deals with the parametrization of § 3 in terms of the Euler angles, which due to spherical 
symmetry are more suitable than the cartesian coordinates. 
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Consider ip, ip, 9 to be the Euler's angles and let 

(p + tfj tp-ip 
«=— . £=— ■ 



The sphere § J can be parametrized as 



tp + ijj 



x\ = cos cos - = cos a cos - , 

2 2 

9 . 9 

x 2 = sin J cos - = sin a cos - , 

2 2 > 





2 




+ v 




2 




-V 




2 


v5 





«i = cos — - — sin - = cos 8 sin - , 
2 2 

■ , ■ 11 

y 2 = sm — - — sin - = sin 8 sin - , 



with < < 7r, — 7T < a < 7r. In the following we shall write the restriction of the one-form 

uj = xidyi - yxdxi + x 2 dy 2 - y 2 dx 2 
to S 3 using Euler's angles. Since 

dxi = — sin a cos ^-da — ^ cos a sin %d9 
2 2 2 

9 1 9 
dx 2 = cos a cos — da — — sin a sin — d9 

9 1 9 
dyi = — sin /3 sin -c?/3 + - cos /3 cos -d9 

9 1 9 

dy 2 = cos /3 sin -c?/3 + - sin /3 cos -d9, 

we obtain 

lo = xidyi — y\dxx + x 2 dy 2 — y 2 dx 2 

( R ■ ■ 6 6 ■ R ■ 9 6 \ r l 

= cos o sm a sm — cos sm p cos a sm — cos — da 

V 2 2 2 2/ 

/ \ 

+ ( — cos a cos - sin (3 sin - + sin a cos - sin - cos (3 ) df3 

\ Z Z Z Z / 

1 / 9 9 9 9 9\ 

+ —( cos a cos - cos 8 cos - + cos a cos 8 sin 2 - + sin a sin 8 cos 2 - + sin 8 sin a sin 2 - ) d9 
2\ 2 2 2 2 2/ 

9 9 9 9 1 

= sin - cos - sin(a — 8)da + sin - cos - sin(a — 8)d8 + — cos(/3 — a)d9 

= - sin 9 sin (da + d 8) + - cos ipd9 



= i ^ sin 9 sin dip + cos ^ g?#^ . 
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Hence 

W|S3 = - ^ sin 9 sin ^ dtp + cos ip d6^j . 

The constraint sin 9 sin f/W^ + cos tpdO — is nonholonomic since 

2du) — cos sin ipd9 /\dp + sin 6* cos ipdip hdip — sin ip dip A d9 ^ . 

We finally arrived at the following characterization of the horizontal curves using Euler's 
angles. 

Lemma 4.1 Let c(s) = (ip(s),ip(s),9(s)) be a curve on § 3 . The curve c is horizontal if and 
only if u)(c) = 0, i.e. 

sm0smipLp + cosip9 = 0. (4.7) 

Before getting into the proof of the connectivity theorem, which is the main result of this 
section, a couple of lemmas are needed. The first one is a standard Calculus exercise. 

Lemma 4.2 Given the numbers a, (3, 7 € M, there is a smooth function 
f : [0, 1] -> M such that 

/(0) = 0, /(l) = a, /'(0)=/3, /'(1)= 7 - 
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Proof: Applying Lemma l4~2l there is a function / : [0, 1] — > M such that 

/(0)=0, /(1)=I, /'(0)=<fo, f'{l) = qi . 

Let q{s) = f'(s). Then f(s) = f* q(u) du and hence / = /(l) = J* q{u) du and q(0) = /'(0) = 
q and q(l) = /'(l) = q x . ■ 

The next result is the theorem of connectivity announced at the beginning of this section. 

Theorem 4.4 Given two points P,Q £ S 3 , there is a smooth horizontal curve joining P and Q. 

Proof: Let P(<fio, V'o, $o) an d Qi'Pi, "4>i> #i) be the coordinates in the Eulerian angles. We shall 
find a horizontal curve c(s) = (ip(s),ip(s), 9(s)) with c(0) = P and c(l) = Q. This is equivalent 
with finding functions ip(s),ip(s),9(s) such that 

sin 0(s) sin ip(s)ip(s) + cos ip(s)0(s) = 0, (4.8) 

which satisfy the boundary conditions 

<K0)=W), #))=V*> 0(0) -0o, 

p(l) = ¥>i, ^(l)=Vi, 0(1) =9,. 
Assuming cos?/>(s) ^ 0, (14. 8|) can be written as 

sin 0(s) tan V>(s)(^(s) + 0(s) = 0. (4.9) 

Let 

9(s) = arcsinp(s), ?p(s) — arctang(s), ip(s) = tpo + s(ipi — tpo) 
for some functions p(s), q(s) which will be determined later. Let k = ipi — tpo- Then (|4.9j) yields 

p'(s) 



p(s)q(s) k 



Separating and solving for p(s) yields 



dp 



pyi -p 2 



-kq(s) ds 



= 0. 



(4.10) 



and integrating, 



tanh 1 



1 



1 

\/l-P 2 (s) 



= tanh 



q(u) du + C\ 
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The constant C\ and the integral q(u) du can be determined from the boundary conditions. 
We obtain 



n/i-p 2 (o) 



tanh(fcCi 



C\ = — tanh" 1 
k 



Vi-p 2 (o) 



v/l-P 2 (l) 



= tanh 



fc( / g(w)du + Ci^) 



f 1 1 

y q(u) du = - 



tanh 



— tanh 1 



v/i-p 2 (i) a/i-p 2 (0) 

where p(Q) = sin 0(0) = sin(9 and p(l) — sin#i. Hence q(s) has to satisfy 



q(u) du 



1 



tanh 1 



1 



COS Ol 



tanh 1 



I cos O | 



The boundary conditions also yield 

g(0)=tanV'O7 q(l) = t&nipi. 

Lemma T4.3I provides the existence of a function q : [0, 1] 
satisfied. 

The function p(s) will be defined by the relation 
1 



(4.11) 



(4.12) 



such that (|4"TTj) and P~T2]) are 



— tanh 



k / q(u) du + tanh 1 



1 



cos 6*o I 



Hence the curve c(s) = ((p(s),ip(s),6(s)) = (tpo + ks, arctang(s), arcsinp(s)) is the desired 
horizontal curve joining the points P and Q on S 3 . ■ 

It is known that in the case of the Heisenberg group if (xq,0) and (xi,0) are two points in 
the {t = 0} plane, then the segment joining the points is a horizontal curve joining the points 
which lies in the {t — 0} plane. The following theorem is an analog for the sphere S 3 . The 
plane will be replaced by the sphere § 2 . 

Theorem 4.5 Given (tpo, tpx) and (9o, 0\), there is a horizontal curve with ip = constant, which 
joins the points with Euler coordinates (<po, ip, 9q) and (ipi,il>,6i). 

Proof: Let (ip, tp, 9) be a horizontal curve with tp = constant. Then equation (|4. T[) can be 
written as 

, dip 1 
-tanV>— = 

dd sm 
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Integrating, yields 

-tanV f^dS = [ -^de^ 
Je d 9 J e sm0 

-taxiip^p(6)-(p(do)J = In tan- - In tan y. (4.13) 

For 9 = 6*i we obtain 

/ tan — \ 

ip = arctan I In /(cpo — pi) I . (4-14) 

V tan — / 

-j 

Solving for p, (|4.13|) yields 

/ 9 9 \ 
— (^o = — In ^ tan — / tan y J / tan ip, 

with -0 denned by (|4.14p . Hence the horizontal curve joining (po, ip, 9o) and (pi,ip, 9\) is 

/ / 9 9 \ \ 

(ip,ip,0) = (ipo - In ( tan-/ tanyj/tan^, 0, 6>J . 

■ 

Let S0 O = {(<y9, -0, #) 6§ 3 ;^ = V'o} be the subset of S 3 with the same Euler angle ip. This 
corresponds to a 2-dimensional sphere. 

Corollary 4.6 Given the points ((po,ipo,6o) and (<fi,ipi,0i) on § 3 , with ipo = ipi, there is a 
horizontal curve in Sij, which joins the points. 



5 The subRiemannian geodesies 

Given two points on the sphere § 3 , we already know that there are horizontal curves joining 
them. The horizontal curve with the shortest length is called a geodesic. Standard arguments 
show that an equivalent problem of finding the geodesies is to find the horizontal curves with 
minimum energy. This means to minimize the action integral 

£l[a 2 (s)+b 2 (s)}ds 

subject to the horizontality constraint (x, y) = (x,y). This will be minimizers given by the 
solutions of the Euler-Lagrange system with Lagrangian 

L(x,x,y,y) = i(a 2 + b 2 ) + X(s)(x 1 y 1 - y x x\ + x 2 y 2 - V2X 2 ), 

where a and b depend on x,x,y,y, see formulas (|3.4[> - (|3.5|> . The function A(s) is a Lagrange 
multiplier function. 
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5.1 The case A(s) = 0. 

The curve which minimizes the energy in this case satisfies the Euler-Lagrange system with the 
Lagrangian L(x, y, x, y) = -a 2 + —b 2 . We have 

dL , dL 

— - = ax 2 + by 2 , q — = -ax 2 - by 2 
ax\ oxi 

dL ,81.,. 

— — = -axi-byi, - — = axi+byi 
ox 2 ax 2 

dL dL 

— - = -ay 2 + bx 2 , —— = ay 2 - bx 2 

oyi dy! 

dL dL . 

—— = ayi-bxi, ——=-ayi + bxi. 

dy 2 dy 2 



Then the Euler-Lagrange system 







d dL 


dL 


d dL 


dL 








ds di\ 


dxi ' 


ds dx 2 


dx 2 ' 








d dL 


dL 


d dL 


dL 








ds dyi 


dyi 


ds dij 2 


dy 2 




becomes 
















ax 2 + by 2 = 


~2(ax 2 




ax\ + by\ = — 2(axx - 


\-byi), 




ay 2 - bx 2 = 


-2(ay 2 - 


- bx 2 ), 


ay i -bxx = -2(ayi - 


- b±i). 


Multiplying the 


first equation by x 2 


, the second by x\ 


, the third by 


y 2 and 


j/i, adding yields 
















a(x\ + x 2 




2 x 2 + yxxx 


- x 2 y 2 - x\y\ 


) 




-2{a(x 2 x 2 


+ X\X\ - 


V V2V2 H 


-yiyx) + b{ 


x 2 y 2 + xtyi - 


y2i 2 - : 




-2b{y 2 x 2 - 


v y\x\ - 


x 2 y 2 - 


xiVi) 







(5.15) 
(5.16) 



a = 0, 

where we have used the horizontality condition (see Proposition 13. 2[) 

2/2*2 + yiii - x 2 y 2 - xxyi = 0, 

the constraint 

x\ + x\ + y\ + y\ = 1 

and its derivative 

x\X\ + x 2 x 2 + y\y\ + y 2 y 2 = 0. 

In a similar way, multiplying the first equation by y 2 , the second by j/i, the third by — x 2 , 
and the fourth by —X\, adding, yields b = 0. We arrived at the following result, which provides 
four constants of motion (three angles and the energy). 
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Theorem 5.1 Given two points P,Q e § 3 , let 7 be a length minimizing curve among all the 
horizontal curves which join P and Q. Then 

1) |-y(s) j is constant along the curve. 

2) The angles between the velocity 7 and the vector fields X , Y and T are constant along the 

curve. 



Proof: 

1) From a = b = it follows that a(s) =constant and b(s) ^constant and hence |7(s)| 

\J a(s) 2 + b(s) 2 = constant. 

2) Using 1), we have 



cos 7, X 
cos 7, Y 



(i,X) = a(s) 

111 \X\ ^a(s) 2 +b(s) 2 

( 7 ,y) b(s) 



i\\Y\ ^/a(s) 2 + b(s) 

tTt = tt/2. 



constant 



= constant 



cost7t= feS = 



l7l \T\ 




Figure 1 : The angle between the velocity 7(s) and the vectors X and Y is constant. 



Remark 5.2 This case when the Lagrange multiplier vanishes yields in the case of the Heisen- 
berg group lines parallel to the x-space. These are particular cases of subRiemannian geodesies. 
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5.2 The case A(s) ^ 0. 
The Lagrangian in this case is 

L = ^a 2 + + \{xiy x + ±22/2 - 2/1^1 - 2/2«2)- 
The Euler-Lagrange system f|5 . 1 5[) - (|5.16[) becomes 

ax 2 + by 2 = -2(a±2 + by 2 + Xyx) - Xyi 

axi+byi = —2{ax\ + byx - Xy 2 ) + Xy 2 

ay 2 -bx 2 = -2(ay 2 - bx 2 + X±x) - Xxx 

ayi — bx\ = —2{ayx — bx\ — \i 2 ) + Xx 2 . 

Multiplying the first equation by x 2 , the second by xx, the third by y 2 and the forth by y\, 
adding yields 

a(x\ + x\ + y\ + y\) + b(y 2 x 2 + y x xx - x 2 y 2 - xxVx) 
=i =o 

= -2(a(x 2 x 2 + xx±i + y 2 y 2 + yxyx) + b(x 2 y 2 + xxyx ~ yix 2 - yxxx) 
V „ ' y ' 

=0 =0 

+X(x 2 yx - xiy 2 + ±xy 2 - i 2 J/i) ) + X{-x 2 yx + x x y 2 - x\y 2 + x 2 y\) <==> 

<■ w ' / " „ ' 

=b =0 

a = -2Xb. 

In a similar way, multiplying the first equation by y 2 , the second by yi, the third by —x 2 
and the forth by —x\, adding we obtain 

a(x 2 y 2 + xxyx - x 2 y 2 - xxyx) + b{y\ + y\ + x\ + x\) 

" v ' » ' 

=0 =1 

= -2{a(x 2 y 2 + ±xy\ - y 2 x 2 - yxXx) + b(y 2 y 2 + yxyx + x 2 x 2 + XxXx) 

=0 

+X{yxV2 - hVi ~ xxx 2 + x 2 Xx) + X(-yiy 2 + y 2 yx - xxx 2 - x 2 xi) 

S v ' " v ' 

=-a =0 

b = 2Xa. 



Hence a and b satisfy the following system 



a = -2Xb (5.17) 
b = 2Xa. (5.18) 
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Multiplying the first equation by a and the second by b and add, yields ad + bb = 0, i.e., 
a 2 + b 2 = r 2 constant along the geodesies. This means that the energy is preserved and the 
velocity of the geodesies have constant length. Let 6{s) be a function such that 

a{s) = r cos 9(s), b(s) — rsm6(s). 

Substituting in the equations (|5.17|) . (|5.18p . we obtain 

sin 0(s)6(s) = 2\{s)sm9(s) 
cos 9(s)9(s) = 2\{s)cos6(s), 

which after dividing the first equation by sin 9 and the second by cos 9 yields 

9{s) = 2A(s). 

The solution is 9(s) = 2A(s) + 9 , where A(s) = / X(u) du. Hence 

Jo 

a(s) = rcos(2A(s) + 6» ), (5.19) 

b(s) = rsm(2A(s)+6» ). (5.20) 

Since A(0) = 0, then a(0) = rcos^o and 6(0) = rsin^o, which provides 9q = tan _1 (6(0)/a(0)). 
Let 7(s) be a subRiemannian geodesic. Since {AT, Y} are orthonormal, 

a = (j,X) = |7| • |*| cos(^X) =rcos(^X) 

=r =1 

and 

b = ( 7 , Y) = | 7 | • \Y\ cos(TTr) = r cos( 7 , Y) = r sin (| - 7 , Y). 

=r =1 

Comparing with (|5.19p . (|5.20[) we arrive at the following result. 

Proposition 5.3 The angles under which the subRiemannian geodesies intersect the integral 
curves of the vector fields X and Y are given by the formulas 

j(s),X j(s) = 2A(s) + 6» , 7 (s),F 7 ( s) = ^ - i(s),X l(s) . 

The main goal now is to find A(s). In order to do this we shall construct an equivalent 
variational problem. We need the following result, which writes the energy in a more friendly 
way. 

Proposition 5.4 If (a;i(s), aegfs), yi(s), J/2(s)j is a subRiemannian geodesic then we have 
± 2 {s) + ± 2 {s) + y 2 (s) + y 2 {s) = a 2 {s) + b 2 {s). 



1G 



Proof: Using the definitions of a and b, the horizontality condition, and the holonomic con- 
straint, we can write the following system 

x\x 2 + 2/12/2 - X2X1 - 2/11/2 = a(s) 

x\yi + x 2 yi - ±22/1 - hxi = b(s) 

±12/1 + ±22/2 - x 2 y 2 - xiyi = 

iixi + X2X2 + 2/12/1 + 2/22/2 = 0. 

This can be written in a matrix way as 



/ x 2 


-X\ 


-2/2 






(xA 




2/2 


-2/1 


x 2 


~X\ 




x 2 


b 


yi 


J/2 


~X\ 


-x 2 




2/i 





\Xl 


x 2 


yi 


2/2 ) 




W 


w 



Since det M = 1, M 1 = M* , M is an orthogonal matrix, and hence preserves the Euclidean 
length of vectors. It follows that 

■ 2 1 • 2 1 '2 1 '2 2,;2 

+ + 2/r + 2/2 = a + & 
along any subRiemannian geodesic, which proves the statement. ■ 

Remark 5.5 The above result holds because the Euclidean metric on K 4 restricted on span{X, Y} 
is the subRiemannain metric. 

The next result will be useful later in the sequel, when we shall find the Lagrange multiplier. 

Proposition 5.6 Let j(s) be a subRiemannain geodesic. Then (7(s),T 7 ( s )) = 7 for any s, 
i.e., the component of the acceleration along the missing direction T — [X,Y] vanishes. 

Proof: Let j(s) — (xi(s),x 2 (s),yi(s),y 2 (s)). Differentiating in the horizontality condition 

xiyi + x 2 y 2 - yixi - y 2 x 2 = 

yields X\y\ + x 2 y 2 — y\x\ — y 2 x 2 = 0, which can be written as 

((x 1 ,x 2 ,y\,y 2 ),(y 1 ,y 2 ,-x 1 ,-x 2 )) = 0, 

or (7, T) =0, where 

T = yid Xl + y 2 d X2 - x x d yi - x 2 d V2 . 
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We shall consider another variational problem, where we replace the energy ~^( a2 + b 2 ) by 
-(i^ + x 2 + yf + y 2 ) an d consider the new Lagrangian 

L*(x,x,y,y) = i(irf +xl + y 2 +A(s)(±ij/i + £22/2 - 2/1 an -mx 2 ). 

The Euler-Lagrange equations provided by the Lagrangian L* describe the same subRieman- 
nian geodesies as the equations associated with the initial Lagrangian L. The Euler-Lagrange 
equations for L* are 

xi = 2Xyi + Xyi 
'±2 = 2Xy 2 + Xy 2 
y\ = — 2\x\ — Xxi 
y 2 = -2A± 2 - Xx 2 . 

Multiplying the first equation by y\, the second by y 2 , the third by —x±, and the fourth by 
— x 2 , adding yields 

2/1 xi + 2/2*2 - xiyi - x 2 y 2 = 2X(y 1 y 1 + y 2 y 2 + x\X\ + x 2 x 2 ) 

v , ' 

=0 

+kyl + yt + xl + xl), 

=1 

which can be written as 

A(s) = 2/1 xx + y 2 x 2 - xxy\ - x 2 y 2 

= ((xi,x 2 , 271,2/2), (2/1,2/2, -a?i, -£2)) 
= (7,T>- 

Using Proposition 15.61 we obtain the following result. 
Proposition 5.7 The Lagrange multiplier A(s) is constant along any subRiemannian geodesies. 
If let A = c/2, with c constant, then Proposition 15.31 yields 

7(s),X 7W = cs + 6» . 

This leads to the following characterization of the subRiemannian geodesies. 

Proposition 5.8 The unit speed curve j(s) is a subRiemannian geodesic on (§ 3 ,AT, Y) if and 
only if the angle between its velocity and the direction of the vector field X increases linearly 
in s. 
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5.3 The group SU(2) 
The special unitary group 

SU(2) = {xi + x 2 i + yj + y 2 k; if + x 2 + y\ + y\ = 1} 
is also a Lie group, which can be identified with the sphere S 3 by the isomorphism 
ip : SU(2) — > § 3 , ip(xtl + x 2 i + j/ij + y 2 k) = (x 1 ,x 2l yi,y 2 ). 

Since 

X q = x 2 d Xl - xid X2 - y 2 d yi + yid V2 

= {x 2 , -xi, -y 2 , yi) ■ V 9 = (x 2 - x x i - y 2 j + j/ x k) • V g 

= -(xjl + x 2 i + yij + y 2 k)i ■ V 9 = -(xi, x 2 , yi, 2/2)1 • V 9 

= -<(<?i),V ? >, (5.21) 



^"g = y 2 d Xl - yid X2 + x 2 d yi - x 1 d V2 

= (V2,-yi,x 2 , -Xi) ■ W q = (y 2 I - yii + x 2 j - xxk) ■ W q 

= -(xil + x 2 i + yij + y 2 k)k ■ V 9 = -(x 1 ,x 2l y 1 ,y 2 )k ■ X7 q 

= -<(gk),V 9 >, (5.22) 

T q = yid Xl + y 2 d X2 - x\d m - x 2 d V2 

= (2/1,2/2, -X!,-x 2 ) ■ V g = {yil + y 2 i - Xij - x 2 k) ■ V q 

= -(xj + x 2 i + yi] + y 2 k)j • V g = ~(x 1 ,x 2 ,y 1 ,y 2 )j ■ V g 

= -((d), V 9 >, (5.23) 

where q = (x%, x 2 , yi, 1/2), V g = (d Xl ,d X2 ,d yi ,d V2 ) and i, j, k are defined in (12. 1|) . It follows 
that the left invariant vector fields X, Y, T, which span the Lie algebra of S 3 correspond on 
SU(2) to — i, k, — j, respectively. 

Remark. An alternative approach to calculate geodesies is to use the Hamiltonian method. 
This method is very complicated in the case of studying geodesies on § 3 , but some information 
we can obtain from the superficial analysis. Using the notations (|5.21ll . (|5.22| . (|5.23|) . the vector 
fields X, Y and T can be written in the form 

X=(-qh,V q ), Y=(-qI 3 ,V 9 ), T = (-qI 2 ,V q ) 

where (•, •) is the usual scalar product in R 4 . Define the Hamilton function 

H = \{X^ + y 2 ) = \ ((qh, 2 + (qh, 2 ) • 
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Then the Hamilton system is the following 

<i = ^r <i = (qiuO-(qh) + (qh,0-(qh) 

f)H (5 - 24) 

i = —g- =*■ i = (Qh,d-(th) + {qh,d-(t;h). 

A geodesic is the projection of the solution of Hamilton's system to the g-space. Since (ql\ , q) = 
{ql2,q} = (qh,q) = 0, multiplying the first equation of (|5.24j) by q we get 

(Si?) = =^ W\ 2 = const. 

We conclude that the solution of the Hamiltonian system belongs to a sphere. Taking the 
constant equals 1 we get geodesies belonging to § 3 . Multiplying the first equation of (]5 . 24[) by 
ql 2 , we get 

(q,qh)=0, (5.25) 

by the role of multiplication between I\ , I 2 , and I3 . The reader easily recognize the horizontality 
condition (x, y) — (x,y) in (|5.25[) . It means that the solution of the Hamiltonian system is a 
horizontal curve. Multiplying the first equation of (|5.24[) by ql\ or by ql 3 , we get 

We see that the Hamilton function can also be written in the form 

H = \((qh,0 2 + (qh,0 2 ) =l((qh,q) 2 + (ql 3 ,q) 2 
If we multiply the first equation of (|5.24[) by q then we get 

\q\ 2 = (qh,0 2 + (qh,0 2 = (qh,q) 2 + (qhA? = 2H. (5.26) 

.12 

Thus the Hamiltonian function give the kinetic energy H = 4j- and it is a constant along 
the geodesies. Notice that (qh,q) = a and {qly,,q) = b by (|3.4[) and (|3.5p . We conclude 
I <7 1 2 = a 2 + b 2 from (|5.26p that corresponds to Proposition This is only the beginning of a 
long term project. There are many problems remain open. For example, given any two points 
on § 3 , how many geodesies connecting them? Is there any abnormal minimizcr in this case? 
What is the action function? What is the volume element which is the solution of a transport 
equation? We will answer these questions in a forthcoming paper. 
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